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Abstract. Let S C {—1,0,1}^ \ {(0,0)}. We address the enumeration of plane lattice 
walks with steps in S, that start from (0,0) and always remain in the first quadrant ■ 
* > 0,j > 0}. A priori, there are 2* problems of this type, but some are trivial. Some 
others are equivalent to a model of walks confined to a half-plane: such models can be solved 
systematically using the kernel method, which leads to algebraic generating functions. We 
focus on the remaining cases, and show that there are 79 inherently different problems to 
study. 

To each of them, we associate a group G of birational transformations. We show that this 
group is finite (of order at most 8) in 2.3 cases, and infinite in the 56 other cases. 

We present a unified way of solving 22 of the 23 models associated with a finite group. 
For all of them, the generating function is found to be D-finite. The 23rd model, known as 
Gessel's walks, has recently been proved by Bostan et al. to have an algebraic (and hence 
D-finite) solution. We conjecture that the remaining 56 models, associated with an infinite 
group, have a non-D-finite generating function. 

Our approach allows us to recover and refine some known results, and also to obtain new 
results. For instance, we prove that walks with N, E, W, SS, SW and NE steps have an 
algebraic generating function. 



1. Introduction 

The enumeration of lattice walks is a classical topic in combinatorics. Many combinatorial 
objects (trees, maps, permutations, lattice polygons, Young tableaux, queues...) can be encoded 
by lattice walks, so that lattice path enumeration has many applications. Given a lattice, 
for instance the hypercubic lattice U^, and a finite set of steps 5 C Z'^, a typical problem is 
to determine how many n-step walks with steps taken from 5, starting from the origin, are 
confined to a certain region A of the space. If A is the whole space, then the length generating 
function of these walks is a simple rational series. If ^ is a half-space, bounded by a rational 
hyperplane, the associated generating function is an algebraic series. Instances of the latter 
problem have been studied in many articles since at least the end of the 19th century [Tl[6]. It 
is now understood that the kernel method provides a systematic solution to all such problems, 
which are, in essence, one-dimensional [IH[2]- Other generic approaches to half-space problems 
are provided in [2Ql [25] . 

A natural next class of problems is the enumeration of walks constrained to lie in the intersec- 
tion of two rational half-spaces — typically, in the quarter plane. Thus far, a number of instances 
have been solved, but no unified approach has emerged yet, and the problem is far from being 
completely understood. The generating functions that have been found demonstrate a more 
complicated structure than those of half-space problems. Some examples are algebraic, but for 
reasons that are poorly understood combinatorially \2Q \ [T3 t [40 ^ [9]. Some examples are, more 
generally, D-finite, meaning that the generating function satisfies a linear differential equation 
with polynomial coefficients |3 H fTS j fTT j I47j. Some examples are not D-finite, having infinitely 
many singularities in the complex plane [TSlBT]. 

We focus in this paper on walks in the plane confined to the first quadrant. The four examples 
of Figure [U illustrate the complexity of this problem: 
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• Kreweras' walks (steps W, SS, and NE): these were first counted in 1965 by Krew- 
eras [37]. He obtained a complicated expression for the number of n-step walks ending 
at (i, j), which simplifies drastically when j = 0. It was then proved by Gessel that the 
associated 3-variable generating function is algebraic [26]. Since then, simpler deriva- 
tions of this series have been obtained [121 US], including an automated proof [36], and 
a purely bijective one for walks ending at the origin [1]. See also [231 [211 [l3] , where the 
stationary distribution of a related Markov chain in the quarter plane is obtained, and 
found to be algebraic. 

• Gessel's walks (steps E, W, NE and SW): Around 2001, Gessel conjectured a sim- 
ple hypergeometric formula for the number of n-step walks ending at the origin. This 
conjecture was proved recently by Kauers, Koutschan and Zeilberger [35]. Even more 
recently, Bostan and Kauers proved that the associated 3-variable generating function 
(which counts all quarter plane walks by the length and the coordinates of the endpoint) 
is in fact algebraic [9]. Strangely enough, the simple numbers conjectured by Gessel had 
not been recognized as the coefficients of an algebraic series before. Both approaches 
involve, among other tools, heavy computer algebra calculations. 

• Gouyou-Beauchamps's walks (steps E, W, NW and SE): Gouyou-Beauchamps 
discovered in 1986 a simple hypergeometric formula for walks ending on the a;-axis |29) . 
We derive in this paper similar expressions for the total number of walks, and for those 
ending at a prescribed position. The associated series are D-finite, but transcendental. 
These walks are related to Young tableaux of height at most 4 [3^ . An afiine deformation 
transforms them into square lattice walks (with N, SS, E and W steps) confined to the 
wedge < j < I. An enumeration of these walks involving the number of visits to the 
diagonal appears in [32l HH |43] . 

• A non-D-finite case is provided by walks with N E, NW and SE steps. Mishna and Rech- 
nitzer estabHshed a compHcated expression for the generating function of these walks, 
from which they were able to prove that this series has infinitely many singularities, and 
thus cannot be D-finite [41] . 

Observe that we have defined 5, the set of steps, as a subset of 1? but that we often use a more 
intuitive terminology, referring to (1,1) as a NE step, for instance. We occasionally abuse the 
coordinate notation and directly write our steps using x's and y's, writing for example xy for a 
SE step. 




Figure 1. Four models of walks in the quarter plane: Kreweras' walks, Gessel's 
walks, Gouyou-Beauchamps's walks, and the non-D-finite example of Mishna 
and Rechnitzer. The numbers qiji) count walks of length n confined to the 
quarter plane that start and end at the origin. We have used the notation 
(a)„ = a(a -I- 1) • • • (a -I- n - 1). 
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Ideally, we seek generic results or combinatorial conditions which assure D-finite (or even 
algebraic) generating functions. One criterion of this type states that if the set S is invariant by 
reflection around a vertical axis (we say for short that it has a vertical symmetry) and consists 
of steps («, j) such that |z| < 1, then the generating function of quarter plane walks with steps in 
S is D-finite [TSl [Tl] . This is also true when the set of steps is left invariant by a Weyl group and 
the walks are confined to a corresponding Weyl chamber [28]. We are not aware of any other 
such criteria. 

1.1. Results 

We restrict our attention to the quarter plane and small steps {i.e., the step set 5 is a subset 
of {—1,0, 1}^ \ {(0,0)}). This includes the four examples above. We first narrow down the 2* 
possible cases to 79 distinct non-trivial problems (Section [2|). Their step sets S are Hsted in 
Tables [l] to [4] in Section [HI These problems fall into two categories, depending on whether a 
certain group associated with S is finite or infinite (Section [3]) . The 23 models associated with 
a finite group turn out to be those that satisfy at least one of the following properties: 

- the step set possesses a vertical symmetry, 

- the vector sum of the vectors in the step set is 0. 

In Section [4] we develop certain general tools that apply to all models — in particular, we 
explain how to write for each of them a functional equation that defines the generating function 
of the walks. Then, we describe a uniform way to solve this equation for all of the models 
associated with a finite group, except one (Gessel's walks). The solutions are all D-finite, and 
even algebraic in three cases (Sections [5] and [6]) . Note that Gessel's walks are also known to 
have an algebraic generating function [9]. We conjecture that all the solutions to models with 
an infinite group are non-D-finite. We conclude in Section [7] with some comments and questions. 
The tables of Section [8] list the 79 models, classified according to the order of the corresponding 
group, and provide references to both the existing literature and to the relevant result of this 
paper. 

1.2. Comments and detailed outline of the paper 

The following technical and/or bibliographical comments may be of interest to readers who 
have already worked on similar problems. One will also find here a more detailed description of 
the contents of the paper. 

The starting point of our approach is a functional equation defining the generating function 
Q{x,y]t) that counts quarter plane walks by the length (variable t) and the coordinates of the 
endpoint (variables x and y) . This equation merely refiects a step by step construction of quarter 
plane walks. For instance, the equation obtained for Kreweras' walks (first example in Figure [T|) 
reads: 

(1 - t{l/x + 1/y + xy)) Qix, y;t)^l- t/x Q(0, y; t) - t/yQ{x, 0; t). 
Note that there is no obvious way to derive from the above identity an equation for, say, (5(0, 0; t) 
or Q{l,l;t). Following Zeilberger's terminology [56], we say that the variables x and y are 
catalytic. One of our objectives is to provide some general principles that may be applied to any 
such Hnear equation with two catalytic variables. The case of linear equations with one catalytic 
variable is well-understood, and the solutions are always algebraic |14) . 

One key tool in our approach is a certain group G{S) of birational transformations that 
leaves the kernel of the functional equation (that is, the coeflicient of Q{x,y;t)) unchanged 
(Section [31) . We have borrowed this group from the little yellow hook by Fayolle, lasnogorodski 
and Malyshev [H], in which the authors study the stationary distributions of Markov chains 
with small steps in the quarter plane. Ever since it was imported from probability theory to 
combinatorics, this group has proved useful (in several disguises, like the obstinate, algebraic, or 
iterated kernel method) to solve various enumeration problems, dealing with walks [TTj [T3j [HI 
[SHI [ini [H] , but also with other objects, like permutations [l^ or set partitions [TTt I55j — the 
common feature of all these problems being that they boil down to solving a linear equation 
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with two catalytic variables. A striking observation, which applies to all solutions obtained so 
far, is that the solution is D-finite if and only if the group is finite. We find that exactly 23 out 
of our 79 quarter plane models give rise to a finite group. 

We then focus on these 23 models. For each of them, we derive in Section [4] an identity 
between various specializations of Q{x, y; t) which we call the orbit sum (or half-orbit sum, when 
there is an x/y symmetry in G{S)). 

In Section m we show how to derive Q{x,y;t) from the orbit sum in 19 out of the 20 models 
that have a finite group and no x/y symmetry. The number of n-step walks ending at 
is obtained by extracting the coefficient of x'^yH'^ in a rational series which is easily obtained 
from the step set and the group. This implies that the generating function Q{x, y; t) is D-finite. 
The form of the solution is reminiscent of a formula obtained by Gessel and Zeilberger for the 
enumeration of walks confined to a Weyl chamber, when the set of steps is invariant under the 
associated Weyl group [28]. Indeed, when the quarter plane problem happens to be a Weyl 
chamber problem, our method can be seen as an algebraic version of the refiection principle 
(which is the basis of [28]). However, its range of applications seems to be more general. We 
work out in details three cases: walks with N, W and SE steps (equivalent to Young tableaux 
with at most 3 rows), walks with N, SS, E, W, SE and NW steps (which do not seem to have 
been solved before, but behave very much like the former case), and finally walks with E, W, 
NW and SE steps (studied in [29]), for which we obtain new explicit results. 

The results of Section [6] may be considered more surprising: For the 3 models that have a 
finite group and an x/y symmetry, we derive the series Q{x,y;t) from the half-orbit sum and 
find, remarkably, that Q{x,y;t) is always algebraic. We work out in details all cases: walks 
with SS, W and NE steps (Kreweras' walks), walks with E, N and SW steps (the reverse steps of 
Kreweras' steps) and walks with N, SS, E, W, NE and SW steps, which, to our knowledge, have 
never been studied before. In particular, we find that the series <5(l,l;t) that counts walks of 
the latter type, regardless of their endpoint, satisfies a simple quartic equation. 

1.3. Preliminaries and notation 

Let A be a commutative ring and x an indeterminate. We denote by A[x] (resp. ^[[x]]) the 
ring of polynomials (resp. formal power series) in x with coefficients in A. If A is a field, then 
A{x) denotes the field of rational functions in x, and A{{x)) the field of Laurent series in x. 
These notations are generalized to polynomials, fractions and series in several indeterminates. 
We denote x = 1/x, so that A[x, x] is the ring of Laurent polynomials in x with coefficients in A. 
The coefficient of x" in a Laurent series F{x) is denoted [x"']F{x). The valuation of a Laurent 
series F{x) is the smallest d such that x'' occurs in F{x) with a non-zero coefficient. 

The main family of series that we use is that of power series in t with coefficients in A[a;,a;], 
that is, series of the form 

n>0,ieZ 

where for all n, almost all coefficients f(i; n) are zero. The positive part of F(x; t) in x is the 
following series, which has coefficients in a;Q[x]: 

[x>]F{x-t) -.^ ^ f{i-n)xH-^. 

n>OA>0 

We define similarly the negative, non-negative and non-positive parts of F{x; t) in x, which we 
denote respectively by [x'^]F{x;t), [x-]F{x;t) and [x-]F{x;t). 

In our generating functions, the indeterminate t keeps track of the length of the walks. We 
record the coordinates of the endpoints with the variables x and y. In order to simplify the 
notation, we often omit the dependence of our series in t, writing for instance Q{x, y) instead of 
Q{x,y]t) for the generating function of quarter plane walks. 

Recall that a power series F{xi, . . . , Xk) € . . . , Xfc]], where K is a field, is algebraic (over 

K(xi, . . . , Xk)) if it satisfies a non-trivial polynomial equation P(xi, . . . , Xfc, i^(a;i, . . . , Xk)) = 0. 
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It is transcendental if it is not algebraic. It is D-finite (or holonomic) if the vector space over 
K{xi, . . . , Xk) spanned by all partial derivatives of F{xi, . . . , Xk) has finite dimension. This 
means that for all z < k, the series F satisfies a (non-trivial) linear differential equation in Xi 
with coefficients in K[xi, . . . , a;^]. We refer to |38[ [39] for a study of these series. All algebraic 
series are D-finite. In Section [5] we use the following result. 

Proposition 1. If F{x,y;t) is a rational power series in t, with coefficients in C{x)[y,y], then 
[y^]F(x,y;t) is algebraic over C{x,y,t). If the latter series has coefficients in C[x,x,y], its 
positive part in x, that is, the series [x-^][y-^]F{x,y;t) , is a 3-variable D-finite series (in x, y 
and t). 

The first statement is a simple adaptation of [25l Thm. 6.1]. The key tool is to expand F(x, y; t) 
in partial fractions of y. The second statement relies on the fact that the diagonal of a D-finite 
series is D-finite (3^. One first observes that there is some k such that [y^]F{x,y;tx'^) has 
polynomial coefficients in x and y, and then applies Remark (4), page 377 of [38] . 

Below we also use the fact that a series F[t) with real coefficients such that [t'^]F{t) ~ Kfj,"n~'^ 
with fc e {1, 2, 3, . . .} cannot be algebraic |22) . 

2. The number of non-equivalent non-simple models 

Since we restrict ourselves to walks with "small" steps (sometimes called steps with small 
variations), there are only a finite number of cases to study, namely 2*, the number of sets 
S formed of small steps. However, some of these models are trivial (for instance 5 = 0, or 
S = {x}). More generally, it sometimes happens that one of the two constraints imposed by 
the quarter plane holds automatically, at least when the other constraint is satisfied. Such 
models are equivalent to problems of walks confined to a half-space: their generating function 
is always algebraic and can be derived automatically using the kernel method [HI [2] . We show 
in Section [2?T] that, out of the 2^ = 256 models, only 138 are truly 2-constraint problems and 
are thus worth considering in greater detail. Then, some of the remaining problems coincide up 
to an x/y symmetry and are thus equivalent. As shown in Section [2?2t one finally obtains 79 
inherently different, truly 2-constraint problems. 

2.1. Easy algebraic cases 

Let us say that a step is a;-positive if z > 0. We define similarly x-negative, y-positive 
and y-negative steps. There are a number of reasons that may make the enumeration of quarter 
plane walks with steps in 5 a simple problem: 

(1) If S contains no x-positive step, we can ignore its i-negative steps, which will never be 
used in a quarter plane walk: we are thus back to counting walks with vertical steps on 
a (vertical) half-line. The solution of this problem is always algebraic, and even rational 
ff 5 = or 5 = {y} or 5 = {y}; 

(2) Symmetrically, if S contains no y-positive step, the problem is simple with an algebraic 
solution; 

(3) If S contains no x-negative step, all walks with steps in S that start from (0, 0) lie in 
the half-plane i > 0. Thus any walk lying weakly above the x-axis is automatically a 
quarter plane walk, and the problem boils down to counting walks confined to the upper 
half-plane: the corresponding generating function is always algebraic; 

(4) Symmetrically, if S contains no y-negative step, the problem is simple with an algebraic 
solution. 

We can thus restrict our attention to sets S containing a;-positive, x-negative, y-positive and 
y-negative steps. An inclusion-exclusion argument shows that the number of such sets is 161. 
More precisely, the polynomial that counts them by cardinality is 

Pi{z) = (1 + ~ 4(1 + 2)5 + 2(1 + zf + 4(1 + zf - 4(H- z) -f 1 

= 2 -h 20 2^ -(- 50 -f- 52 z^ + 28 8 z'^ z^ 
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In the expression of Pi{z), one of the 4 terms (1 + z)^ counts sets with no x-positive step, one 
term (1 + z)^ those with no a;-positive nor x-negative step, one term (1 + z)^ those with no x- 
positive nor y-positive step, and so on. All the sets S we have discarded correspond to problems 
that either are trivial or can be solved automatically using the kernel method. 

Among the remaining 161 sets S, some do not contain any step with both coordinates non- 
negative: in this case the only quarter plane walk is the empty walk. These sets are subsets of 
{x, y, xy, xy, xy}. But, as we have assumed at this stage that S contains a;-positive and y-positive 
steps, both xy and xy must belong to S. Hence we exclude 2^ of our 161 step sets, which leaves 
us with 153 sets, the generating polynomial of which is 

P^{z) = Pi{z) - z^{l + zf = + 17 z^ + 47 z" + 51 z^ + 28 z^ + 8 z^ + z^. 

Another, slightly less obvious, source of simplicity of the model is when one of the quarter plane 
constraints impHes the other. Assume that all walks with steps in S that end at a non-negative 
abscissa automatically end at a non-negative ordinate (we say, for short, that the x-condition 
forces the y-condition) . This implies in particular that the steps y and xy do not belong to S. 
As we have assumed that S contains a y-negative step, xy must be in S. But then x cannot 
belong to S, otherwise some walks with a non-negative final abscissa would have a negative 
final ordinate, like x followed by xy. We are left with sets S C {x,y,xy,xy,xy} containing 
xy, and also xy (because we need at least one x-positive step). Observe that these five steps 
are those lying above the first diagonal. Conversely, it is easy to realize that for any such set, 
the x-condition forces the y-condition. The generating polynomial of such super- diagonal sets 
is z^(l -I- z)^. Symmetrically, we need not consider sub-diagonal sets. An inclusion-exclusion 
argument reduces the generating polynomial of non-simple cases to 

Psiz) = P2{z) - 2z'^{l + zf + z'^ = llz34-41z'' + 49z^ + 28z^ + 82^ + z^, 

that is to say, to 138 sets S. 

2.2. Symmetries 

The eight symmetries of the square act on the step sets. However, only the x/y symmetry 
(refiection across the first diagonal) leaves the quarter plane fixed. Thus two step sets obtained 
from one another by applying this symmetry lead to equivalent counting problems. As we want 
to count non-equivalent problems, we need to determine how many among the 138 sets S that 
are left have the x/y symmetry. We repeat the arguments of the previous subsection, counting 
only symmetric models. We successively obtain 

- = (1 + z)2(l + ;^2)3 _ 2(1 + z)(l + z2) + 1, 

P^y-" = Pf - z2 = 3z3 -t- 5z4 + 5z5 + 4z6 + 2z^ + z^. 

For instance, the term we subtract from P^^™ to obtain Pj^™ counts symmetric subsets of 
{x, y, xy, xy, xy} containing xy and xy. The generating polynomial of (inherently different) 
models that are neither trivial, nor equivalent to a 1-constraint problem is thus 

i (P3 + PT'^) = 7 z^ + 23 z* + 27 z'^ 16 z^ 5 z^ z*^. 
This gives a total of 79 models, shown in Tables [T] to HI 

3. The group of the walk 

Let 5 be a set of small steps containing x-positive, x-negative, y-positive and y-negative steps. 
This includes the 79 sets we wish to study. Let S{x, y) denote the generating polynomial of the 
steps of S: 

S{x,y)= xW- (1) 
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It is a Laurent polynomial in x and y. Recall that x stands for 1/x, and y for 1/y. Let us write 

S{x,y) = A.i{x)y + Ao{x) + Ai{x)y = B^i{y)x + Bo{y) + Bi{y)x. (2) 

By assumption, Ai, Bi, A^i and B^i are non-zero. Clearly, S{x,y) is left unchanged by the 
following rational transformations: 

^ , , ( B^Ay) \ , , , ^ / -A.Ax) 
$ : (x, y) X y and V : (x, y) i-^ x, y 



Si (2/) "7 M{x) 

Note that both $ and ^' are involutions, and thus foirational transformations. By composition, 
they generate a group that we denote G(5), or G if there is no risk of confusion. This group 
is isomorphic to a dihedral group D„ of order 2n, with n e N U {oo}. For each g G G, one has 
S{g{x,y)) = S{x,y). The sign of 5 is 1 (resp. —1) if g is the product of an even (resp. odd) 
number of generators $ and 

Examples 

1. Assume S is left unchanged by a reflection across a vertical line. This is equivalent to saying 
that S{x,y) = S{x,y), or that Bi{y) = i?_i(y), or that Ai{x) = Ai{x) for i = —1,0,1. Then 
the orbit of {x,y) under the action of G reads 

(x, y)^{x, y)^{x,C{x)y)'r^{x,C{x)y)'r^{x, y), 

with C{x) = so that G is flnite of order 4. 

Note that there may exist rational transformations on (x, y) that leave S{x, y) unchanged but 
are not in G. For instance, if 5 = {N, , E, W}, the map (x, y) 1-^ (y, x) leaves S{x, y) unchanged, 
but the orbit of (x, y) under G is {(x, y), (x, y), (x, y), (x, y)}. 

2. Consider the case S = {x,y,xy}. We have A-i{x) = x, Ai{x) — 1, i3_i(y) — 1, i3i(y) = y. 
The transformations are 

$ : (x, y) (xy, y) and * : (x, y) (x, xy), 
and they generate a group of order 6: 

(x, y)^(xy, y)^(xy, x)^(y, x)^(y, xy)^(x, xy)^(x, y). (3) 

3. Consider now the case S = {x, y, xy}, which differs from the previous one by a rotation of 90 
degrees. We have A„i(x) = 1, Ai{x) ~ x, _B_i(y) = 1, Bi{y) = y. The two transformations are 

$ : (x, y) (xy, y) and * : (x, y) (x, xy), 

and they also generate a group of order 6: 

(x, y)^(xy, y)^(xy, x)^(y, x)^(y, xy)^(x, xy)^(x, y). 

As shown by the following lemma, this is not a coincidence. 

Lemma 2. Let S and S be two sets of steps differing by one of the 8 symmetries of the square. 
Then the groups G{S) and G{S) are isomorphic. 

Proof. The group of symmetries of the square is generated by the two reflections A (across 
the first diagonal) and V (across a vertical Hue). Hence it suffices to prove the lemma when 
S = A(5) and when S — V{S). We denote by $ and ^ the transformations associated with S, 
and by $ and ^' those associated with S. 

Assume S ~ A(S'). We have Ai{x) = Bi[x) and i?i(y) = Ai{y). Denote by 5 the involution 
that swaps the coordinates of a pair: 5[x,y) = (y,x). An elementary calculation gives 

$ = (5 o vj/ o (5 and = 5 o o 5 

so that the groups G{S) and G{S) are conjugate by 5. 
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Assume now S — V{S). We have Ai{x) ~ Ai{x) and Bi{y) = B^i{y). Denote by v the 
involution that replaces the first coordinate of a pair by its reciprocal: v{x,y) = {x,y). An 
elementary calculation gives 

<(> = w o $ o w and ~ v o o v 
so that the groups G{S) and G{S) are conjugate by v. ■ 

Theorem 3. Out of the 79 models under consideration, exactly 23 are associated with a finite 
group: 

- 16 have a vertical symmetry and thus a group of order 4, 

- 5 have a group of order 6, 

- 2 have a group of order 8. 

Proof. Given that $ and 4" are involutions, the group they generate is finite of order 2n if and 
only if := 5" o $ has finite order n. It is thus easy to prove that one of the groups G{S) has 
order 2n: one computes the rrfi^ iterate 0™ for 1 < m < n, and checks that only the last of 
these transformations is the identity. 

We have already seen in the examples above that models with a vertical symmetry have a 
group of order 4. We leave it to the reader to check that the models of Tables [2] and [3] have 
groups of order 6 and 8, respectively. These tables give the orbit of {x, y) under the action of G, 
the elements being listed in the following order: (x, y), $(a;, y), ^ o ^{x, y), and so on. 

Proving that one of the groups G{S) is infinite is a more difficult task. We apply two different 
strategies, depending on S. The first one uses valuations and works for the five step sets of 
Figure [2j These are the sets of our collection for which all elements (i, j) satisfy i + j > 0. We 
are very grateful to Jason Bell, who suggested to us a second strategy which turned out to apply 
to the remaining cases. 



Figure 2. Five step sets with an infinite group. 
1. The valuation argument 

Let z be an indeterminate, and let x and y be Laurent series in z with coefficients in Q, of 
respective valuations a and b. We assume that the trailing coefficients of these series, namely 
[z°]x and [z^]y, are positive. Let us define x' by ^{x,y) = {x',y). Then the trailing coefficient 
of x' (and y) is positive, and the valuation of x' (and y) only depends on a and b: 

( (-a + b{v'i"l-v^{'^),b) if&>0, 
(j)(a,b) := (va\(x'),va\(y)) = < > , , , ^ ( 

^ ' ^ ^' ^" \ {-a + b{d^^\-d[y^),b) if6<0, 

where v'f^ (resp. d^^^) denotes the valuation (resp. degree) in y of Bi(y), for i = ±1. Similarly, 
^'(a;, y) := (x, y') is well-defined, and the valuations of x and y' only depend on a and b: 



V'(a,6) := (val(x),val(y')) = < /' ,^ ,(x) 



a,-b + a{v^_^l-v[''^)) ifa>0, 
fe + a(rfL1 -4""^)) ifa<0. 



where w'^"* (resp. d'f'^) denotes the valuation (resp. degree) in x of Ai{x), for i — ±1. 

In order to prove that G is infinite, it suffices to prove that the group G' generated by (j) and 
■0 is infinite. To prove the latter statement, it suffices to exhibit (a, b) G such that the orbit 
of (a, b) under the action of G' is infinite. 
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Let S be one of the five sets of Figure [21 Then A^i{x) ~ x and B^i{y) = y, so that 
w^i = rf*^! = I'^i = rf^i = 1- Also, v[^^ = v[^^ = — 1 as 5 contains the steps xy and xy. Hence 
the transformations (j> and ip read: 

[ {-a + 2b,b) if6>0, [ {a,2a-b) if a > 0, 



It is easy to check, by induction on n > 0, that 

(V'o0)"(l,2) = (2n+l,2n + 2) and (^(V' o 2) = (2?i + 3, 2?! + 2). 

(All these pairs have positive entries, so that we never need to know dj^'' or Hence the 

orbit of (1,2) under the action of (p and ip is infinite, and so are the groups G' and G. 

We believe, from our computer experiments, that the groups G' associated with the remaining 
51 models of Table [Hare finite, and hence, cannot be used to prove that G is infinite. Instead, 
we use for these models a different argument based on the fixed points of 9 = ^' o $. 

2. The fixed point argument 

We are left with 51 models. Thanks to Lemma [2l we only need to prove that (roughly) a quarter 
of them are associated with a finite group: if G'(5) is infinite, then G{S) is infinite for all sets 
S that differ from 5 by a symmetry of the squarqj. This leaves 14 models to study, listed in 
Table [51 

Assume = 5* o $ is well-defined in the neighborhood of (a, 6) e C^, and that this point is 
fixed by O. Note that a and b are algebraic over Q. Let us write Q — (0i, 82), where Qi and 
62 are the two coordinates of 6. Each 9i sends the pair {x,y) to a rational function of x and 
y. The local expansion of 6 around (a, b) reads 

e{a + u,b + v) = (a, b) + (u, v)JaM + 0{u'^) + O(w^) + 0{uv), 

where Ja,b is the Jacobian matrix of 8 at (a, b): 

( aei 962, \ 



Ja,b — 

Iterating the above expansion gives, for m > 1, 

6'"(a + u,b + v) = (a, b) + {u, v)J^^ + 0{u^) + 0{v'') + 0{uv). 

Assume G{S) is finite of order 2n, so that 6 has order n. Then Q"'{a + u,b + v) = (a, b) + (u, v), 
and the above equation shows that J"^, is the identity matrix. In particular, all eigenvalues of 
Ja,b are roots of unity. 

This gives us a strategy for proving that a group G(S) is infinite: find a fixed point (a, 5) 
for 6, and compute the characteristic polynomial xi^) of the Jacobian matrix Ja,b- This is a 
polynomial in X with coefficients in Q{a,b). In order to decide whether the roots of x are roots 
of unity, we eliminate a and b (which are algebraic numbers) from the equation x(X) = to 
obtain a polynomial x{^) ^ Qi^] that vanishes for all eigenvalues of Ja,b- if none of its factors is 
cyclotomic, we can conclude that G{S) is infinite. As all cyclotomic polynomials of given degree 
are known, this procedure is effective. 

Let us treat one case in detail, say S = {x,y,y,xy} (the first case in Table [5]). We have 

6(a;, y) o $(a;, y) = {xy, x + y) . 



Why a quarter, rather than an eighth? Recall that, if two (distinct) models differ by an x/y symmetry, only 
one of them appears in Tabled] 
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Every pair (a, b) such that a'^ + — 1 and b = Xjc? is fixed by 9. Let us choose one such pair. 
The Jacobian matrix reads 




and its characteristic polynomial is 

X{X) ~ dct(X \A-Ja.h) = + X(l + a*) + + . 
Let X be a root of this polynomial. By eliminating a (which satisfies a** + = 1), we obtain 

X(X) := + 9 + 31 + 62 X-^ + 77 X* + 62 + 31 + 9 X + 1 = 0. 

This polynomial is irreducible, and distinct from all cyclotomic polynomials of degree 8. Hence 
none of its roots are roots of unity, no power of Ja,fc is equal to the identity matrix, and the 
group G{S) is infinite. 

This strategy turns out to work for all models of TablelH This table gives, for each model, the 
algebraic equations defining the fixed point (a, b) that we choose (for instance, the "condition" 
a^ + a'' — 1 occurring on the first line means that a^ + a^ — 1 = 0), and a polynomial x(X) g Q[-'^] 
that vanishes at all eigenvalues of the Jacobian matrix, in factored form. One then checks that 
no factor of this polynomial is cyclotomic. 

It may be worth noting that this second strategy does not work for the five models of FigureO 
in the first three cases, 9 has no fixed point; in the last two cases, it has a fixed point (a, 6), 
but the sixth power of the Jacobian matrix Ja,h is the identity (of course, 6^ is not the identity; 
more precisely, the expansion of 9^(a + m, 6 + v) involves cubic terms in u and v). g 

Remark. To put this discussion in a larger framework, let us mention that the group of 
birational transformations or (or of the projective plane P^(C)), called the Cremona group, has 
been the object of many studies in algebraic geometry since the end of the 19th century |34[I54). 
It seems possible that, given the attention already paid to the classification of finite subgroups of 
this group (see e.g. [HE]), there is a generic or automatic test for finiteness that can be applied 
to all our examples. 



4. General tools 

Let S be one of the 79 step sets of Tables[T]to|ll Let Q be the set of walks that start from (0, 0), 
take their steps from S and always remain in the first quadrant. Let q{i,j;n) be the number 
of such walks that have length n and end at position (i, j). Denote by Q{x,y;t) = Q{x,y) the 
associated generating function: 

Q{x,y;t)= q{i,j;n)x'yH''. 
It is a formal power series in t with coefficients in j/]. 

4.1. A FUNCTIONAL EQUATION 

Lemma 4. As a power series in t, the generating function Q{x, y) = Q{x, y; t) of walks with 
steps taken from S starting from (0, 0) and staying in the first quadrant is characterized by the 
following functional equation: 

K{x, y)xyQ{x, y) = xy - txA^i{x)Q{x, 0) - tyB.i{y)Q% y) + ieQ(0, 0) 

where 

K{x, y) = 1 - tS{x, y)^l-t ^ x'y^ 

(ij)es 

is called the kernel of the equation, the polynomials A^i{x) and B^i{y) are the coefficients of y 
and X in S{x,y), as described by and e is \ if (—1, —1) is one of the allowed steps, and 
otherwise. 
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Proof. We construct walks step by step, starting from the empty walk and concatenating a 
new step at the end of the walk at each stage. The empty walk has weight 1. The generating 
function of walks obtained by adding a step of S at the end of a walk of Q is tS{x,y)Q{x,y). 
However, some of these walks exit the quadrant: those obtained by concatenating a y-negative 
step to a walk ending at ordinate 0, and those obtained by concatenating an x-negative step 
to a walk ending at abscissa 0. Walks ending at ordinate (resp. abscissa) are counted by 
the series Q{x,Q) (resp. (5(0, y)). Hence we must subtract the series tyA-i{x)Q{x,Q) and 
txB-i(y)Q(0,y). However, if (—1,-1) € S, we have subtracted twice the series counting walks 
obtained by concatenating this step to a walk ending at (0,0): we must thus add the series 
etxyQ {0,0). This inclusion-exclusion argument gives 

Q{x, y) = l + tS{x, y)Q{x, y) - tyA^^{x)Q{x, 0) - txB_^{y)Q{0, y) + etxyQ{0, 0), 

which, multiplied by xy, gives the equation of the lemma. 

The fact that it characterizes Q{x,y;t) completely (as a power series in t) comes from the 
fact that the coefficient of t" in Q(x, y; t) can be computed inductively using this equation. This 
is of course closely related to the fact that we have used a recursive description of walks in Q to 
obtain the equation. B 

4.2. Orbit sums 

We have seen in Section[3]that all transformations g of the group G associated with the step set 
S leave the polynomial S{x, y) unchanged. Hence they also leave the kernel K{x, y) = l — tS{x, y) 
unchanged. Write the equation of Lemma [4] as 

K{x, y)xyQ{x, y) ^ xy ~ F{x) - G{y) + teQ{0, 0), 

with F[x) = txA_i{x)Q{x,0) and G(y) = tyB.i{y)Q{0,y). Replacing {x,y) hy<^>{x,y) = {x',y) 
gives 

K{x, y)x'yQ{x', y) = x'y ~ F{x') - G{y) + teQ{0, 0). 

The difference between the former and latter identities reads: 

K{x, y) {xyQ{x, y) - x'yQ{x', y)) = xy - x'y - F{x) + F{x'). 

The term G{y) has disappeared. We can repeat this process, and add to this identity the 
functional equation of LemmalU evaluated at {x',y') = 'if{x',y). This gives: 

Kix, y) {xyQix, y) - x'yQ{x', y) + x'y'Q{x', y')) ^ xy ~ x'y + x'y' - F{x) - G{y') + teQ(0, 0). 

Now the term F{x') has disappeared. If G is finite of order 2n, we can repeat the procedure 
until we come back to (4* o (f>)"(a;, y) = {x, y). That is to say, we form the alternating sum of the 
equations over the orbit of {x,y). All unknown functions on the right-hand side finally vanish, 
giving the following proposition, where we use the notation 

forgGG, giAix,y))~Aig{x,y)). 

Proposition 5 (Orbit sums). Assume the group G{S) is finite. Then 

sign(g)5(a;yQ(a;, y; t)) = — - ^ sign{g)g{xy). (4) 
sec ^^^'2^'^^ gee 

Observe that the right-hand side is a rational function in a;, y and t. We show in Section [5] 
that this identity implies immediately that 19 of the 23 models having a finite group have a 
D-finite solution. 

The 4 remaining models are Gessel's model {x, x, xy, xy} (which we do not solve in this paper) 
and the three models with steps Si = {x, y, xy}, S2 — {x, y, xy}, and 5 = 5i U ^2. For each of 
these three models, the orbit of {x, y) is 

{x, y)^{xy, y)^{xy, x)^{y, x)^{y, xy)^{x, xy)^{x, y). 
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and exhibits an xjy symmetry. That is, (y, x) belongs to the orbit of {x,y). Moreover, if 
g{{x,y)) = iy,x), then sign((7) = —1. Thus the right-hand side of ^ vanishes, leaving 

xyQ{x, y) - xQ{xy, y) + yQ{xy, x) = xyQ{y, x) - xQ{y, xy) + yQ{x, xy). 

But this identity directly follows from the obvious relation Q(x, y) = Q(jj, x) and does not bring 
much information. In Section[6l we solve these three obstinate models by summing the functional 
equation over one half of the orbit only. Given that A-i{x) ~ B-i{x), the identity resulting 
from this half-orbit summation reads as follows. 

Proposition 6 (Half-orbit sums). Denote Si = {x,y,xy} and S2 = {x,y,xy}. If the set of 
steps is Si, S2 or 5i U ^2, then 

. . xy-x + y~ 2txA^i{x)Q{x, 0) + teQ(0, 0) 

xyQ{x, y) - xQ{xy, y) + yQ{xy, x) = — . 

K{x,y) 

Remark. For Gessel's walks, the orbit of {x,y) is shown in Tabled Proposition [5] reads 

X! sign{g)g{xyQ{x,y)) = 0, 

sec 

although no obvious symmetry explains this identity. 

4.3. The roots of the kernel 

Recall that the kernel of the main functional equation (Lemma [IJ is 

K{x,y) = l-t ^'y'- 

Lemma 7. Let 

A{x) = (1 - tAo{x)f - 4t^A^i{x)Ai{x). 

Let Yo{x) and Yi{x) denote the two roots of the kernel K{x,y), where K{x,y) is seen as a 
polynomial in y. These roots are Laurent series in t with coefficients in Q(x); 



l-tAo(:r)-7AR 1 - tAo{x) + ^fMx) 

= ' ^^^-"^^ 2^1^) ■ 

Their valuations in t are respectively 1 and —1. Moreover, 1/K{x,y) is a power series in t with 
coefficients in Q[x,x,y,y], and the coefficient of y^ in this series can be easily extracted using 

(6) 



K{x,y) ^/A{^ \l ~ yYoix) l-y/Yi{x) 

Proof. The equation K{x, Y) = also reads 

Y = t (A_i(.t) + YAo{x) + Y^Aiix)) . (7) 

Solving this quadratic provides the expressions of Yo{x) and Yi{x) given above. As A(a;) = 
1 + 0{t), the series Yi has valuation —1 in t, and first term l/{tAi{x)). The equation 

Yo{x)Yiix)^^-^ 
Ai{x) 

then shows that Yo{x) has valuation 1. This is also easily seen on ([7]), which in turn implies that 
Yq has coefficients in Q[a;,a;]. The equation 

1 Aq{x) 



Yoix)+Yiix) 



tAi{x) Ai{x) 



shows that for n > 1, the coefficient of t" in Yi{x) is also a Laurent polynomial in x. This is not 
true, in general, of the coefficients of t~^ and 
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The identity ((6]) results from a partial fraction expansion in y. Note that both Yq and 1/Fi 
have valuation 1 in t. Hence the expansion in y of 1/K{x,y) gives 



Yoix)- 



^y'^K{x,y) ^ 



if j<0 

(8) 

if j>0. 



4.4. Canonical factorization of the discriminant A(a;) 

The kernel can be seen as a polynomial in y. Its discriminant is then a (Laurent) polynomial 
in x: 

A(a;) = (1 - Mo(x))2 - U"^ A^i{x)Ai{x). 

Say A has valuation —5, and degree d in x. Then it admits 5 + d roots Xi = Xi{t), for 
1 < i < d + d, which are Puiseux series in t with complex coefficients. Exactly 6 of them, say 
Xi, . . . , Xg, are finite (and actually vanish) at t = 0. The remaining d roots, ^5+1, . . . , Xg^d, 
have a negative valuation in t and thus diverge at i = 0. (We refer to |52[ Chapter 6] for 
generalities on solutions of algebraic equations with coefficients in C{t).) We write 

A(a;) = AoA_(x)A+(a;), 



with 



s 



A_(i) = A_(x;t) = l[{l-xX,), 

i=l 
S+d 

A+{x) = A+{x;t) = J] (l-x/X,), 



and 



i=6+l 

S+d 



Ao ^ Ao{t) = (-1)^^^ = i-lAx'Mx) n X,- 

lli=l^' i=5+l 

It can be seen that Aq (resp. A_(a;), A+(a;)) is a formal power series in t with constant term 1 
and coefficients in C (resp. C[x], C[a;]). The above factorization is an instance of the canonical 
factorization of series in Q[x, x][[t]], which was introduced by Gessel and has proved useful 
in several walk problems since then \W\ \TE\ [T3] . It will play a crucial role in Section [6l 

5. D-FINITE SOLUTIONS VIA ORBIT SUMS 

In this section, we first show that 19 of the 23 models having a finite group can be solved 
from the corresponding orbit sum. This includes the 16 models having a vertical symmetry, plus 
3 others. Then, we work out the latter 3 models in details, obtaining closed form expressions for 
the number of walks ending at prescribed positions. 

5.1. A GENERAL RESULT 

Proposition 8. For the 23 models associated with a finite group, except from the four cases 
S ~ {x,y,xy}, S = {x,y,xy}, S = {x,y,x,y,xy,xy\ and S = {x,x,xy,xy\, the following 
holds. The rational function 

^^""'^'^^ " K{x,y;t) ^ "^g^(5)5(a^y) 

is a power series in t with coefficients in Q(a;) [y, y] (Laurent polynomials in y, having coefficients 
in Q{x) ). Moreover, the positive part in y of R{x, y; t), denoted R^{x, y; t), is a power series in 
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t with coefficients in Q[x,x,y]. Extracting the positive part in x of R^{x,y;t) gives xyQ{x,y;t) . 
In brief, 

xyQ{x,y;t) = [x>][y>]R{x,y;t). (9) 
In particular, Q(x,y;t) is D-finite. The number of n-step walks ending at [i,]) is 

q{i,3- n) = [x'+'y=+^] sign(ff)5(a;2;) j S{x, y)" (10) 

where 

S{x,y)= ^ x^y'. 
(p,«)e5 

Proof. We begin with the 16 models associated with a group of order 4 (Table [T]). We will then 
address the three remaining cases, S = {x, y, xy}, S = {x, x, xy, xy} and S = {cc, x, y, y, xy, xy}. 

All models with a group of order 4 exhibit a vertical symmetry. That is, K{x,y) = K{x,y). 
As discussed in the examples of Section [3l the orbit of {x, y) reads 

(x, y)^{x, y)*^{x,C{x)y)^{x,C{x)y)'^{x, y), 

with C{x) = ^7^. The orbit sum of Proposition [5] reads 

xyQ{x, y) - xyQ{x, y) + xyC{x)Q{x, C{x)y) + xyC{x)Q{x,C{x)y) = R{x, y). 

Clearly, both sides of this identity are series in t with coefficients in Q(x)[?/, y]. Let us extract 
the positive part in y: only the first two terms of the left-hand side contribute, and we obtain 

xyQ{x, y) - xyQ{x, y) = i?+(x, y). 

It is now clear from the left-hand side that R^{x, y) has coefficients in <Q[x, x, y]. Extracting the 
positive part in x gives the expression ^ for xyQ{x,y), since the second term of the left-hand 
side does not contribute. 

Let us now examine the cases S — {x,y,xy}, S = {x,x,xy,xy} and S = {x,x,y,y,xy,xy}. 
For each of them, the orbit of {x, y) consists of pairs of the form {x°'y'', x'^y'^) for integers a, b, c and 
d (see Tables [2] and [3]) . This implies that R{x,y) is a series in t with coefficients in Q[x,x,y,y]. 
When extracting the positive part in x and y from the orbit sum of Proposition \E\ it is easily 
checked, in each of the three cases, that only the term xyQ{x,y) remains in the left-hand side. 
The expression of xyQ{x,y) follows. 

Proposition [T] then impHes that Q{x,y;t) is D-finite. The expression of q{i,j;n) follows from 
a mere coefficient extraction. B 

5.2. Two MODELS WITH ALGEBRAIC SPECIALIZATIONS: {x,y,xy} AND {x, X, y , y, xy , xy} 

Consider the case S = {x,y,xy} = {W, N,SE}. A walk if with steps taken from S remains 
in the first quadrant if each of its prefixes contains more N steps that SE steps, and more SE 
steps than W steps. These walks are thus in bijection with Young tableaux of height at most 3 
(Figure [3]), or, via the Schensted correspondence [49], with involutions having no decreasing 
subsequence of length 4. The enumeration of Young tableaux is a well-understood topic. In par- 
ticular, the number of tableaux of a given shape — and hence the number of n-step walks ending 
at a prescribed position — can be written in closed form using the hook-length formula [l9]. It is 
also known that the total number of tableaux of size n and height at most 3 is the nth Motzkin 
number |48]. Here, we recover these two results (and refine the latter), using orbit and half-orbit 
sums. Then, we show that the case S = {x,x,y,y,xy,xy}, which, to our knowledge, has never 
been solved, behaves very similarly. In particular, the total number of n-step walks confined to 
the quadrant is also related to Motzkin numbers (Proposition [TO]) . 
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N-N-SE-N-SE-W-W-N-SE-W 



Figure 3. A Young tableau of height 3 and the corresponding quarter plane 
walk with steps in {W, N, SE}. 

Proposition 9. The generating function of walks with steps W, N, SE confined to the quar- 
ter plane is the non-negative part (in x and y) of a rational series in t having coefficients in 
Q[x,x,y,y]: 

Q{x,y;t) = [x^y^]R{x,y;t), with Rix,y;t) ^ ^- \ ,_ ^ -^ ■ 

1 — t(x + y + xy) 

In particular, Q(x, y) is D-finite. The number of walks of length n = 3m + 2i + j ending at 

IS 

. ^ (^ + l)(J^ + l)(^+j + 2)(3m + 2^+j)! 
m!(m + z + l)!(m + z+j + 2)! 

In particular, 

9(0,0; 3m) = ,, ^ V3^, 

m!(m + lj!(m + 2j! irm^ 

so that Q{0,0;t) — and hence Q{x,y;t) — is transcendental. 

However, the specialization Q{x, \/x;t) is algebraic of degree 2; 



^, _ , 1-tx- ^/l~2xt + f^x^ ~ Af^x 
Q[x, 1/x; t) = 



2xt^ 

In particular, the total number of n-step walks confined to the quadrant is the n^^ Motzkin 
number: 



n>0 k=0 



2k V k 



Proof. The orbit of {x, y) under the action of G is shown in |[3]) . The first result of the proposition 
is a direct application of Proposition [HI with R{x,y) — R{x,y)/{xy). It is then an easy task to 
extract the coefficient of x^yH" in R{x,y;t), using 

r iT,- _N„ (3m + 2i + ;■)! 

[xY] ix + y + xyr = ^ .... ... 

m\[m + i)\[m + I + jy. 

if 71 = 3m + 2i + j. 

The algebraicity of Q{x,x) can be proved as follows: let us form the alternating sum of the 
three equations obtained from Lemma [4] by replacing {x,y) by the first 3 elements of the orbit. 
These elements are those in which y occurs with a non-negative exponent. This gives 

K{x,y){xyQ{x,y) - xy^Q{xy,y) + x^yQ{xy,x)) = xy - xy'^ + x'^y - tx^Q{x,{)) - txQ{Q,x). 

We now specialize this equation to two values of y. First, replace y by x: the second and third 
occurrence of Q in the left-hand side cancel out, leaving 

K{x, x)Q{x, x) = l- tx^Q[x, 0) - fxQ(0, x). 

For the second speciaHzation, replace y by the root Yo[x) of the kernel (see This is a 

well-defined substitution, as 10(2;) has valuation 1 in t. The left-hand side vanishes, leaving 

= xYq{x) - xYQ{xf -f x^Yo{x) - tx^Q{x, 0) - txQ{0,x). 
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By combining the last two equations, we obtain 

K{x,x)Q{x,x) = I - xYo{x) + xYo{xf ~ x^Yoix). 
The expression of Q(x,x) follows. 

■ 

The case S = {N, , W, E, SE, NW} is very similar to the previous one. In particular, the orbit 
of [x, y) is the same in both cases. The proof of the previous proposition translates almost 
verbatim. Remarkably, Motzkin numbers occur again. 

Proposition 10. The generating function of walks with steps N, SS, W, E, SE, NW confined to 
the quarter plane is the non-negative part (in x and y) of a rational function: 

r ^ ~ , ~, X (I — xy) {\ — x^y) (l — xy'^) 

Q{x,y-t)^[x^y^]R{x,y-t), with R{x,y;t) = / _ 

1 — t[x + y + x + y + xy + xy) 

In particular, Q(x,y;t) is D-finite. The specialization Q{x, \/x;t) is algebraic of degree 2; 



^, , I . l-tx-tx+J{l- t{x + x)Y - 4<2(1 + x){l + x) 
Q(x,l/x;t) = 



2*2 (1 + a;) (1 + .t) 

:e total number of n-step wt 

Motzkin number. 



In particular, the total number of n-step walks confined to the quadrant is 2" times the n^^ 



g(M;t).i^^^^V^?pIiIM.Vr'^' 2" fn^f2k 



8*2 ^ ^ k+l\2kl\k 

n>0 k=0 \ / \ 

5.3. The case S = {x,x,xy,xy} 

Consider quadrant walks made of E, W, NW and SE steps. These walks are easily seen to be 
in bijection with pairs of non-intersecting prefixes of Dyck paths: To pass from such a pair to 
a quadrant walk, parse the pair of paths from left to right and assign a direction (E, W, NW 
or SE) to each pair of steps, as described in Figure IH Hence the number of walks ending at a 
prescribed position is given by a 2-by-2 Gessel-Viennot determinant [22] and we can expect a 
closed form expression for this number. 




Figure 4. A pair of non- intersecting prefixes of Dyck paths corresponding to 
the quarter plane walk E-W-E-E-NW-SE-W-E-E-NW. 

Also, note that the linear transformation {i + maps these quarter plane walks 

bijectively to walks with E, W, N and SS steps confined to {(«, j) : < j < i}. In this form, 
they were studied by Gouyou-Beauchamps, who proved that the number of n-step walks ending 
on the X-axis is a product of Catalan numbers [l^. His interest in these walks came from a 
bijection he had established between walks ending on the a;-axis and Young tableaux of height 
at most 4 [30]. This bijection is far from being as obvious as the one that relates tableaux of 
height at most 3 to quarter plane walks with W, N and SW steps (Section l5.2p . 

Here, we first specialize Proposition [8] to obtain the number of n-step walks ending at (i, j) 
(Proposition [TT|) . Then, we perform an indefinite summation on i, or j, or both i and j to obtain 
additional closed form expressions, including Gouyou-Beauchamps's (Corollary [12]). 
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Proposition 11. The generating function of walks with steps E, W, NW, SE confined to the 
quarter plane is the non-negative part (in x and y) of a rational function: 

Q{x,y;t) = [x~y~]R{x,y;t), 

with 

~ {l-x)il + x){l-y){l-x^y)il-xy)il + xy) 
R{x,y;t) = — -— — _ . 

1 — t\x + a; + + xy) 

In particular, Q{x, y; t) is D-finite. The number of walks of length n = 2m + i ending at {i, j) is 

,. . ^ (i + l)(j + l)(i + j + 2)(i + 2j + 3) /2TO + i + 3\ /2™ + i + 3^ 
gu, 7; Zra + i) = — -r- 

^ ' (2m + i + l)(2m + i + 2)(2m + i + 3)2 \^ m-j )\ m + 1 

In particular, 

6(2m)!(2m + 2)! 24 ■ 4^- 
9(0,0; 2m) = — — — — — — — — 

m!(m + l)!(m. + 2)!(m + 3)! TTm° 

so that (5(0, 0;t) — and hence Q{x,y;t) — is transcendental. 

Before we prove this proposition, let us perform summations on i and j. Recall that a hy- 
pergeometric sequence (/(fc))fe is Gosper summahle (in k) if there exists another hypergeometric 
sequence {g{k))k such that f{k) = g{k + 1) — g{k). In this case, indefinite summation can be 
performed in closed form [46l Chapter 5]: 

fel 

f{k) - 9{h + 1) - 9{ka). 

k—ko 

The numbers q{i,j; n) of Proposition [TT] have remarkable Gosper properties, from which we now 
derive Gouyou-Beauchamps's result for walks ending on the x-axis, and more. 

Corollary 12. The numbers q{i,j;n) are Gosper-summable in i and in j . Hence sums of the 
formJ2iL^^^q{i,j;n) andJ2'f=j^q{i,j]n) have closed form expressions. In particular, the number 
of walks of length n ending at ordinate j is 

(j + l)(2j7i)!(2m + 2)! 

if n = 2m, 



qi-J;^) := ^5(i, j;ri) = < 

i>0 



(m - j)!(m + l)P{m+j + 2)! 
2(j + l)(2m+l)!(2m + 2)! 



(m - j)!(m + l)!(m, + 2)!(m + j + 2)! 
Similarly, the number of walks of length n — 2m + i ending at abscissa i is 

(i + l)(2m + i)!(2m + i + 2)! 



if n = 2m + 1. 



q{i,-;2m + i) := Yq{i,j;2m + i) = 



m\{m + l)!(m + i + l)!(m + i + 2)!' 

In particular, as many 2m-step walks end on the x- and y-axes: 

(2m)!(2m + 2)! 



g(-,0;2m) =9(0, -;2m) 



?n!(m+ l)!2(7Ti + 2)!' 

The numbers q{i, — ;n) and q{—,j;n) defined above are again Gosper-summable in i and j re- 
spectively. Hence sums of the form 

il il jl jl 

^g(i,-;n) = ^^g(i,j;n) and ^ g(-, j: n) = ^ ^ g(i, j; n), 

i=io i-k) j j=jo j=jo i 
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which count walks ending between certain vertical or horizontal lines, have simple closed form 
expressions. In particular, the total number of quarter plane walks of length n is 



q{-,-;n) := ^ q{i,j;n) = 
t.j>a 



' (2?7i)!(2m + 1)! 
mP{m + 1)P 

(2m+ l)!(2m + 2)! 
TO!(TO + l)!2(m + 2)! 

The asymptotic behaviours of these numbers are found to be 



if n = 2m, 
if n — 2m + 1. 



which shows that the series Q{l,l;t), Q{0,l;t) and Q{l,0;t) are transcendental. 

Proof of Proposition]!^ and Corollary]!^ The orbit of {x, y) under the action of G, of cardinaHty 
8, is shown in Table[3l The expression of Q(x, y; t) given in Proposition [11] is a direct apphcation 
of Proposition [SI with R{x,y) = R{x,y)/{xy). It is then an easy task to extract the coefficient 
of x^yH^^ in R{x, y; t), using 



X + X + xy + xy ^ (1 + y){x + xy) and [x'^y-']{x + x -\- xy + xy)" = 



m + I J \m — J 



for n = 2m + i. This proves Proposition [TTl 

For the first part of corollary, that is, the expressions of q{—,j; n) and q{i, — ; n), it suffices 
to check the following identities, which we have obtained using the implementation of Gosper's 
algorithm found in the sumtools package of Maple: 



with 

92{hj;m) 



q{2i,j;2m) = gi{i, j;m) - gi{i + 1, j;m), 
q{2i + l,j;2m+l) = g2ii,j;m)-g2{i + l,j;m), 
q{i,j;2m + i) = 33(2, j; m) - 53(4, j + 1; m), 

2{l+j){m + l + 2i{i + l+j)){2m)l (2to+1)! 
(m — i — j)! (m — i + 1)\ (m + i + 1)! (m + i+j + 2)1' 
2(l+j)(m + l+j(l + 2z) + 2 (l+i)^)(2m + l)! (2m + 2)! 

{m-i- j)l {m-i + l)\{m + i + 2)\ {m + i+j + 3)1 
{1 + i) {m + 1 + {i + j + 1) {1 + j)) {2m + iy. {2m + i + 2)\ 



{m-j)\ (m + 1)! (m + i + 2)! (m + i+j + 2)! 

These three identities respectively lead to 

9(-,i;2m) =.gi(0,j;m), g(-, j; 2m + 1) = 52(0, j; m) and g(i, - ; 2m + i) = ^3(1, 0; m) 

as stated in the corollary. 

For the second part, we have used the following identities, also discovered (and proved) using 
Maple: 

q{2i,-;2m) = 54(1; m) - 34(1 + 1; m), 

q{2i + l,-;2m+l) = 35(1; m) - 55(1 + 1; m), 

q{-J;2m) = ge{j ; m) - ge{j + 1; m) , 

q{-,j;2m + l) = gr{j]m) ~ gr{j + l;m), 
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with 



96ij;m) 
97{j-:rn) 



(2 to)! (2to + 1)! 



(to — i)! (to — i + 1)! (to + «)! (to + i + 1)! ' 

(2to+ 1)! (2 TO + 2)! 
(to - i)! (to - i + 1)! (to + I + 1)! (to + i + 2)! ' 

(2 to)! (2to + 1)! 

(to - j)!TO! (to + 1)! {171+ j + 1)!' 
(2to+ 1)! (2to + 2)! 



(to-j)! (to + 1)! (to + 2)! (TO + j + 1)!' 

Note that these identities give two ways to determine the total number of n-step walks in the 
quadrant, as 

(/(—j — ; 2m) = 94(0; to) = 36(0; to) and (?(—,— ; 2to + 1) = 95(0; to) = 57(0; to). 



6. Algebraic solutions via half-orbit sums 

In this section we solve in a unified manner the three models whose orbit has an x/y symmetry: 
Si = {x, y, xy}, S2 = {x, y, xy} and 5 = 5i U 52. Remarkably, in all three cases the generating 
function Q(x, y; t) is found to be algebraic. Our approach uses the algebraic kernel method 
introduced by the first author to solve the case S = Si, that is, Kreweras' model |13t Section 2.3]. 
We refer to the introduction for more references on this model. The case S = S2 was solved by 
the second author in [1^, and the case S ~ Si^J S2 is, to our knowledge, new. 

Recall from Proposition [6] that for each of these three models, 

. ^ xy-x + y~ 2txA_i{x)Q{x, 0) + teQjO, 0) 

xyQ{x, y) - xQ[xy, y) + yQ[xy, x) = — r . 

K{x,y) 

Extract from this equation the coefficient of in the left-hand side, only the second term 
contributes, and its contribution is xQd{x), where Qd{x) = Qd{x]t) is the generating function 
of walks ending on the diagonal: 



Qd{x;t) = ^ f^x'q{i,i;n). 



n.i>0 



The coefficient of j/° in the right-hand side can be easily extracted using ^ . This gives 

-xQd{x) = 1— [xYo{x)-x + -^-2txA.i{x)Q{x,0)+teQiO,0) 
or, given the expression ^ of Yq and the fact that YoYi = x, 

xQdix) = f x{l-tAoix)) ___ 2tncA_i{x)Q{x,Q)+teQ{Q,0) 



Let us write the canonical factorization of A(a;) = AoA_|_(x)A_ (a;) (see Section [4~4)) . Multiplying 
the equation by Ai{x)\J /S.-{x) gives 

./Kl¥)(^^-xAi{x)Qd{x)) = 

1 ' •^*^-^^^'^°'^^ll_i^i(a;)-2te^_i(a;)Ai(a;)Q(a;,0)+teAi(a;)g(0,0) ) . (11) 



v/AoA+(x-) V t 

Each term in this equation is a Laurent series in t with coefficients in Q[a;,a;]. Moreover, very 
few positive powers of x occur in the left-hand side, while very few negative powers in x occur 
in the right-hand side. We will extract from this equation the positive and negative parts in a;. 
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and this will give algebraic expressions for the unknown series Qd{x) and Q{x, 0). From now on, 
we consider each model separately. 

6.1. The case S = {x,y,xy} 

We have A-i{x) = 1, Ao{x) = x, Ai{x) = x and e = 0. The discriminant A(a;) reads 
(1 — tx)^ — At^x. The curve A(a;;i) = has a rational parametrization in terms of the series 
W = W{t), defined as the only power series in t satisfying 

W ^t{2 + W^). (12) 

Replacing t by W/{2 + W^) in A{x) gives the canonical factorization as A(a;) = AoA+(x) A_ (x) 
with 



^ ' ^ ' Ax 4a;2 



Ao--^, l^+{x)^l~xW\ A_(a;) = l 

Extracting the positive part in x from ifTTj) gives 

x_ {2t^x'^Q{x,Q)-x + 2t)W W 
1 ^ 2tVl + T' 

from which we obtain an expression of Q{x, 0) in terms of W . Extracting the non-positive part 
in X from (fTTj) gives 



/ WiW^ + A) fx ^ , A X W 

from which we obtain an expression of Qd{x). We recover the results of |13) . 

Proposition 13. Let W = W{t) be the power series in t defined by (fT2| . Then the generating 
function of quarter plane walks formed of\N, SS and NE steps, and ending on the x-axis is 



Qix,0;t) 



tx \2t X \W x) ^ ) 



Consequently, the length generating function of walks ending at (i,0) is 

where Ci = + 1) is the i-th Catalan number. The Lagrange inversion formula gives the 

number of such walks of length m = 3m + 2i as 

/■no n-N 4'"(2i + l) f2i\f3m + 2i 

g(i,0;3m + 2i) - ^ ' ' 



{m + i + l){2m + 2i + l)\i J \ m 

The generating function of walks ending on the diagonal is 

^ , . W-x 
Qd(X] t) = — , + X t. 

ty^l - xW{l + W^/A) + x^W^/4 

Note that Q{0, 0) is algebraic of degree 3, Q{x, 0) is algebraic of degree 6, and Q{x, y) (which 
can be expressed in terms of Q[x,Q), and Q{0,y) = Q{y,0) using the functional equation we 
started from) is algebraic of degree 12. 

6.2. The case S = {x,y,xy} 

The steps of this model are obtained by reversing the steps of the former model. In particular, 
the series Q(0, 0) counting walks that start and end at the origin is the same in both models. 
This observation was used in |40] to solve the latter case. We present here a self-contained 
solution. 
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We have A^i{x) = x, Ao{x) — x, Ai{x) ~ 1 and e — 1. The discriminant A(2;) is now 
(1 — tx)^ — At^x, and is obtained by replacing a; by a; in the discriminant of the previous model. 
In particular, the canonical factors of A(a;) are 

Ao = -^, A+(x-)^l - ^ x+^x\ ^A.x) = l~xW\ 

where W = W{t) is the power series in t defined by fT2|). Extracting the coefficient of in (fTT|) 
gives 

W"^ _ W {W^ + 4 + 8i g(0, 0)) 
2t ^ 16< 

from which we obtain an expression of Q(0, 0). Extracting the non-negative part in x from (fTT|) 
gives 

X _ {2xt^Q{x,Q)~xt^Q{0,0) + t-x^ +x^t)W W 

1: 2t ~ 2xtVl - xW{W^ + 4)/4 + x'^W'^/A 2xt 

from which we obtain an expression of 0) . Finally, extracting the negative part in x from (fTTI) 

gives 

X Qd (x) \ X ^ _ W 



t X J \ X t 2t 2xt^ 

from which we obtain an expression of Qd[x). We have thus recovered, and completed, the 
results of [40]. Note in particular how simple the number of walks of length n ending at a 
diagonal point is. 

Proposition 14. Let W = W{t) be the power series in t defined by (fT2l) . Then the generating 
function of quarter plane walks formed o/ N, E and SW steps and ending on the x-axis is 



W{4:~- W^) t-x'^+ tx^ (2 2;2 - xW'^ - W) ^Jl - xW{W^ + 4)/4 + x'^W'^/A 
Q{x,0;t) = — . 

The generating function of walks ending on the diagonal is 

^ , , xW{x + W)-2 1 

Qd{x]t) = ^ , ' - H ^. 

^ ^ 2ta;Vl - xW^ tx^ 

Consequently, the length generating function of walks ending at is 

w2i+l /o-\ 

[^']Qd{x; t) = -TTTT-^ . (2^ + 4 - (2z + l)W^) . 



4*+it(i + 2) \ i 

The Lagrange inversion formula gives the number of such walks of length n = 3m + 2i as 

4™(i + l)2 f2i + l\f?,m + 2i^ 

q[i, i; 3m + 2i) — 



(m + j + l)(2m + 2i + l)\ i J \ m 

Note that Q(0, 0) is algebraic of degree 3, Q{x, 0) is algebraic of degree 6, and Q{x, y) (which 
can be expressed in terms of Q{x, 0), Q(0, y) = Q{y, 0) and Q{0, 0) using the functional equation 
we started from) is algebraic of degree 12. 

6.3. The CASB S ~ {x,y,x,y,xy,xy} 

We have A^i{x) = 1+5;, Ao{x) = x + x, Ai{x) ^ I + x and e = 1. The discriminant A(a;) 
is now (1 — t{x + i))^ — 4t^(l + x){l + x), and is symmetric in x and x. Two of its roots, say 
Xi and X2, have valuation 1 in t, and the other two roots are 1/Xi and I/X2. By studying the 
two elementary symmetric functions of Xi and X2 (which are the coefficients of A_(a;)), we are 
led to introduce the power series Z = Z{t), satisfying 

1 - 2 Z + 6 Z2 - 2 ^3 + ^4 

O^z)^ (1^) 
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and having no constant term. Replacing t by its expression in terms of Z in A(x) provides the 
canonical factors of A(x) as 

Ao = ^, A+ix) = l-2Z-^±^x + Z'x', A_(x) = A+(x). 

As in the previous case, extracting from (fTT|) the coefficient of gives and expression of Q(0, 0): 

Q(0 0) = 

Extracting then the positive and negative parts of ifTTj) in x gives expressions of Q{x, 0) and 
Qd{x). 

Proposition 15. Let Z = Z{t) he the power series with no constant term satisfying l|13p . and 
denote ^ 

A+(x) = l-2Z^i±|^x + ZV. 

Then the generating function of quarter plane walks formed of N, SS, E, W, SE and NW steps, 
and endinq on the x-axis is 



{Zjl -Z) + 2xZ -{l-Z) x^) TaTM 
' 2teZ(l-Z)(l + x)2 

Z(l - Zf + Z {Z^ + iZ^ -5Z + 2) X- {1-2Z + 7Z^ - AZ^) x^ + x^Z{l - Zf 

2txZ{l - Z)2(l + x)2 
The generating function of walks ending on the diagonal is 

I- Z -2xZ + x'^Z{Z -I) 1 



tx{\ + x){Z^l)^/K^ tx{l+x) 



Note that Q(0, 0) is algebraic of degree 4, Q(x, 0) is algebraic of degree 8, and Q(x, y) (which 
can be expressed in terms of Q{x, 0), Q(0, y) = 0) and Q(0, 0) using the functional equation 
we started from) is algebraic of degree 16. However, Q = Q(l, 1) has degree 4 only, and the 
algebraic equation it satisfies has a remarkable form: 

Q{l + tQ)il + 2tQ + 2t^Q'^) = 

1 — DC 

Also, Motzkin numbers seem to be lurking around, as in Proposition [TOl 
Corollary 16. Let N = N{t) he the only power series in t satisfying 

N = t{l + 2N + 'iN^). 

Up to a factor t, this series is the generating function of the numhers 2"M„, where Mn is the 
Motzkin number: 

L"/2J on 



^>0 fc=0 

Then the generating function of all walks in the quadrant with steps N, S, E, W, SE and NW is 



and the generating function of walks in the quadrant ending at the origin is 

0(0,0;,) = 

The proof is elementary once the algebraic equations satisfied by (3(1,1) and Q(0,0) are 
obtained. 
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7. Final comments and questions 

The above results raise, in our opinion, numerous natural questions. Here are some of them. 
The first two families of questions are of a purely combinatorial, or even bijective, nature ("explain 
why some results are so simple"). Others are more closely related to the method used in this 
paper. We also raise a question of an algorithmic nature. 



7.1. Explain closed form expressions 

We have obtained remarkable hypergeometric expressions for the number of walks in many 
cases (Propositions [9] to [HI) . Are there direct combinatorial explanations? Let us underHne a 
few examples that we consider worth investigating. 

— Kreweras' walks and their reverse: The number of Kreweras' walks ending at (i, 0) 
is remarkably simple (Proposition [T3|) . A combinatorial explanation has been found 
when i = 0, in connection with the enumeration of planar triangulations [4]. To our 
knowledge, the generic case remains open. If we consider instead the reverse collection 
of steps (Proposition[T3]), then it is the number of walks ending at (i, i) that is remarkably 
simple. This is a new result, which we would like to see explained in a more combinatorial 
manner. 

— Motzkin numbers: this famous sequence of numbers arises in the solution of the cases 
iS = {x,y,xy} and S = {x,x,y,y,xy,xy} (Propositions [9l and [T0| . The first problem is 
equivalent to the enumeration of involutions with no decreasing subsequence of length 4, 
and the occurrence of Motzkin numbers follows from restricting a bijection of Frangon 
and Viennot [21] ■ The solution to the second problem is, to our knowledge, new, and 
deserves a more combinatorial solution. Can one find a direct explanation for why the 
respective counting sequences for the total number of walks of these two models differ 
by a power of 2? Is there a connection with the 2" phenomenon of |19j? 

— Gessel's walks: although we have not solved this case {S = {x,x,xy,xy}) in this 
paper, we cannot resist advertising Gessel's former conjecture, which has now become 
Kauers-Koutschan-Zeilberger's theorem: 

,(0,0;2n).16"(^/^)"(l/2)" 



(5/3)„(2), 



Certain other closed form expressions obtained in this paper are less mysterious. As discussed 
at the beginning of Section [531 walks with E, W, NW and SE steps are in bijection with pairs 
of non-intersecting walks. The Gessel- Viennot method (which, as our approach, is an inclusion- 
exclusion argument) expresses the number of walks ending at as a 2-by-2 determinant, 
thus justifying the closed form expressions of Proposition [TTJ The extension of this theory by 
Stembridge [53] allows one to let i, or j, or both i and j vary, and the number of walks is now 
expressed as a pfaffian. Hence the closed forms of Corollary [TH are not unexpected. However, 
one may try to find direct proofs not involving the inclusion-exclusion principle. Moreover, the 
following question may be interesting per se: 

— Walks with E, W, NW and SE steps: can one explain bijectively why as many 2m-step 
walks end on the x- and y-axes? Recall that those ending on the x-axis were counted 
bijectively in [29j . 

Another well-understood case is that of quarter plane walks with N, E, SS and W steps. The 
number of such walks ending at the origin is a product of Catalan numbers, and this has been 
explained bijectively, first in a recursive manner [18], and more recently directly, using again 
certain planar maps as intermediate objects |5]. Another argument, based on the refiection 
principle and thus involving minus signs, appears in [ST] and applies to more general endpoints. 
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7.2. Explain algebraic series 

A related problem is to explain combinatorially, via a direct construction, why the three 
models of Section[6]have algebraic generating functions. Given the connection between Kreweras' 
walks and planar triangulations [1], this could be of the same complexity as proving directly that 
families of planar maps have an algebraic generating function. (Much progress has been made 
recently on this problem by Schaeffer, Di Francesco and their co-authors.) And what about 
Gessel's walks (with steps E, W, NE and SW), which we have not solved in this paper, but have 
very recently been proved to have an algebraic generating function as well [9]? 

7.3. Models with a vertical symmetry 

When S is invariant by a reflection across a vertical axis, the group G{S) has cardinality 4 and 
Proposition [8] gives the generating function Q{x, y; t) as the positive part of a rational function. 
We have not worked out the coefficient extraction in any of these 16 cases. This may be worth 
doing, with the hope of finding closed form expressions in some cases. However, according to [8], 
there is little hope to find an algebraic solution for Q{l,l;t). 

7.4. Models with an infinite group 

Two of the 56 models that are associated with an infinite group (Table \^ have been proved 
to have a non-D-finite generating function |41]. We conjecture that this holds for all models 
with an infinite group. This conjecture is based on our experimental attempts to discover a 
differential equation satisfied by the generating function, and much strengthened by the further 
attempts of Bostan and Kauers (Hj, which are based on the calculation of 1000 terms of each 
generating function. It also relies on the fact that all equations with two catalytic variables and 
an infinite group that have been solved so far have a non-D-finite solution. How could one prove 
this conjecture, for instance in the case S = {N, E, NE, SW}? And what is the importance, if any, 
that can be attributed to the fact that the step sets with finite groups either exhibit a symmetry 
across the y-axis or have a vector sum of zero? 

7.5. Automatic derivation of differential equations 

Our D-finite but transcendental solutions are expressed as the positive part (in x and y) 
of a 3- variable rational series in t, x and y (Proposition [8]). Can one derive automatically 
from these expressions differential equations satisfied by (5(0, 0;t) and (3(1, 1;^)? This would 
be a convenient way to fill in the gap between our work and the paper [8], where differential 
equations are conjectured for the series (5(1, 1; i)- 

For the algebraic solutions of Section [6l it is easy to derive from our results first an algebraic 
equation satisfied by (5(0, 0;t) (or (5(1, l;i)), and then a differential equation satisfied by this 
series, using the Maple package Gfun |50) . 

7.6. Variations and extensions 

It is natural to ask to which similar problems the approach used in this paper could be 
adapted. To make this question more precise, let us underline that such problems may involve, 
for instance, putting weights on the walks, allowing more general steps, or considering higher 
dimensions. However, the very first question is whether Gessel's model can be dealt with using 
the ideas of this paper! 

Weighted paths. One may try to adapt our approach to solve refined enumeration problems. 
For instance, some authors have studied the enumeration of walks in a wedge, keeping track 
not only of the number of steps, but also of the number of contacts, or visits to the boundary 
lines [32 ^ l44 j |43]. Of course, other statistics could be considered. 

Another natural way to add weights, of a more probabilistic nature, consists in studying 
Markov chains confined to the quarter plane. The weight of a walk is then its probability. An 
entire book is devoted to the determination of the stationary distributions of such chains |21) . 
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These distributions are governed by functional equations similar to ours, but without the length 
variable t (since only the stationary regime is considered). This difference makes the problem 
rather different in nature, and indeed, the tools involved in [21] are much more analytic than 
algebraic. A natural way to set the problem back in the algebraic playground (to which our 
power series methods belong) is to keep track of the length of the trajectories, which boils down 
to studying the law of the chain at time n. This was done in [13] for a probabilistic version of 
Kreweras' walks, using a variant of the method presented in this paper. An asymptotic analysis 
of the solution should then yield the limiting/stationary distribution. This was however not done 
in [13]. Instead, we enriched our algebraic approach with a few basic analytic arguments to solve 
directly the equation that describes the stationary distribution. This solution is in our opinion 
more elementary than the original ones [HI [23]. It is worth noting that the bivariate series that 
describes the stationary distribution is algebraic, and that this Kreweras chain is actually the 
main algebraic example of [21]. In view in the results presented in this paper, it is natural to ask 
if one could design probabilistic versions of the other three algebraic models (Propositions [14] 
and [HI plus Gessel's model) that would also yield algebraic stationary distributions. 

Conversely, it is natural to ask whether certain tools from [21], other than the group of the 
walk, could be adapted to our power series context. We are thinking in particular of the material 
of Chapter 4, which is devoted to the case where the group of the walk is finite, and (under this 
hypothesis) to the conditions under which the stationary distribution has an algebraic generating 
function. 

More general steps. The fact that we only allow "small" steps plays a crucial role in our ap- 
proach, and more precisely in the definition of the group of the walk (Section [3]). However, this 
does not mean that models with larger steps are definitely beyond reach. First, it is always pos- 
sible to write a functional equation for the series Q{x, y; t), based on a step-by-step construction 
of the walk. If no step has a coordinate smaller than —1, the right-hand side of the equation only 
involves Q{0,y;t) and Q{x,0;t), but otherwise more unknown functions, depending exclusively 
on X or y, appear. Another important difference with the present setting is that the kernel has 
now degree larger than 2 (in x or y). One can still define a group, but acting on pairs (x, y) than 
cancel the kernel. We refer to [H] for the solution of a simple example, with steps (2,-1) and 
(-1,2). 

Higher dimension. Finally, a natural question is to address 3-dimensional problems as those 
studied experimentally in [8]. Provided one focusses on walks with small steps, the key ingre- 
dients of our approach — the functional equation and the group of the walk — can indeed be 
adapted in a straighforward manner to this higher-dimensional context. 

8. Tables 

The tables below Hst the 79 step sets S we consider, classified according to the cardinality of 
the group G{S). The first three tables contain sets for which G{S) has cardinality 4, 6 and 8 
respectively. The orbit of (.t, y) under the action of this group is listed in the second column. The 
third column lists the steps of S. The fourth one displays the numbers q{l, l;n) and g(0,0;n) 
that respectively count all quarter plane walks and quarter plane walks ending at the origin. 
We have given the reference of these sequences when they appear in the Encyclopedia of Integer 
Sequences [51]. In the rightmost column, we give references on this model, both in this paper 
and in other papers. 
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# 



G{S) 



S 



g(l,l;n) 
9(0, 0;n) 



References 



(x,y),(x, y), 



1, 2, 6, 18, 60, 200, 700, 2450, 8820, 31762 (A005666) 

1, 0, 2, 0, 10, 0, 70, 0, 588, 0, 5544, 0, 56628 
(A005568) 



1, 1, 4, 9, 36, 100, 400, 1225, 4900, 15876 (A018224) 
1, 0, 1, 0, 4, 0, 25, 0, 196, 0, 1764, 0, 17424 (A001246) 



1, 2, 10, 39, 210, 960, 5340, 26250, 148610, 761796 
1, 0, 2, 0, 18, 0, 255, 0, 4522, 0, 91896, 0, 2047452 



1, 3, 18, 105, 684, 4550, 31340, 219555, 1564080 



1, 0, 3, 6, 38, 160, 905, 4830, 28308, 166992 



Prop.ig [m [HI [g EH [31] 



Prop, [a [ni[2Hi[i3 



Prop. [HI [in [28] 



Prop. [HI 



{x,y), {x,y), 



1, 1, 3, 7, 19, 49, 139, 379, 1079, 3011 



0, 0, 0, 2, 0, 0, 0, 28, 0, 0, 0, 



0, 0, 0, 20020 



Prop. [HI 



2, 9, 34, 151, 659, 2999, 13714, 63799, 298397, 
D8415, 6678827 



Prop, [a [U 



0, 1, 3, 4, 20, 65, 175, 742, 2604, 9072, 36960,139392 



{x,y), {x,y), 



2, 7, 23, 84, 301, 1127, 4186, 15891, 60128, 230334 



0, 1, 0, 4, 0, 20, 0, 126, 0, 882, 0, 6732 



3, 15, 74, 392, 2116, 11652, 64967, 365759, 2074574 



1, 0, 2, 3, 12, 40, 145, 560, 2240, 9156, 38724, 166320, 
728508 



Prop, [a [U 



Prop. [HI [U 



10 



{x,y), {x,y), 

y a+i+s)' 
y x+i+x) 



1, 2, 8, 29, 129, 535, 2467, 10844, 50982, 231404 



1, 0, 1, 0, 6, 0, 55, 0, 644, 0, 8694, 0, 128964 



1, 3, 16, 86, 509, 3065, 19088, 120401, 771758 



1, 0, 2, 3, 20, 60, 345, 1400, 7770, 36876, 204876 



Prop. [HI [U 



Prop. [HI [n] 



11 



12 



{x,y), ix,y), 
{x,y{x+l + x)), 
{x,y{x+l + x)) 



1, 1, 3, 6, 17, 34, 121, 265, 969, 2246, 8351, 20118 



1, 0, 1, 0, 4, 0, 20, 0, 126, 0, 882, 0, 6732 



1, 2, 7, 26, 105, 444, 1944, 8728, 39999, 186266 



1, 0, 2, 3, 12, 40, 145, 560, 2240, 9156, 38724, 166320, 
728508 



Prop. [HI [n] 



Prop. [HI [n] 
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13 
14 


{x,y), {x,y), 

\-^iy x+x ) 




X 




1, 1, 5, 13, 61, 199, 939, 3389, 16129, 61601, 295373 


Prop.ig [U 

Prop, [a [U 


1, 0, 1, 0, 6, 0, 55, 0, 644, 0, 8694, 0, 128964 

1, 2, 11, 49, 277, 1479, 8679, 49974, 301169, 1805861 




1, 0, 2, 3, 20, 60, 345, 1400, 7770, 36876, 204876 


15 
16 


{x,y), {x,y), 
{x,y{x + x)), 
{x,y{x + x)) 




/ \ 




1, 1, 2, 3, 8, 15, 39, 77, 216, 459, 1265, 2739, 7842 


Prop, [a [IB So] 
Prop. [HI [n] 


1, 0, 0, 0, 2, 0, 0, 0, 28, 0, 0, 0, 660 

1, 2, 6, 21, 76, 290, 1148, 4627, 19038, 79554, 336112 


1, 0, 1, 3, 4, 20, 65, 175, 742, 2604, 9072, 36960 



Table 1: The group G{S) is isomorphic to -D2- These 16 models 
have a D-finite generating function. 



# 


GiS) 


S 


q{l,l;n) 
g(0,0;n) 


References 


1 

2 


{x,y), {xy,y), 
{xy,x), {y,x), 
{y, yx), (x,yx) 






1, 1, 2, 4, 9, 21, 51, 127, 323, 835, 2188, 5798 
(A001006) 


Prop. El [28] 






1, 0, 0, 1, 0, 0, 5, 0, 0, 42, 0, 0, 462 (A005789) 

1, 2, 8, 32, 144, 672, 3264, 16256, 82688 (A129400) 


Prop.fTOl [28] 


1, 0, 2, 2, 12, 30, 130, 462, 1946, 7980, 34776, 153120 


3 
4 
5 


{x,y), {xy,y), 
{xy,x), {y,x), 
{y,xy), {x,xy) 




n 


1, 1, 3, 7, 17, 47, 125, 333, 939, 2597, 7183 


prop.H [sails] 

Prop. [12 [IQ] 
Prop.[l5| 


1, 0, 0, 2, 0, 0, 16, 0, 0, 192, 0, 0, 2816 (A006335) 
1, 2, 4, 10, 26, 66, 178, 488, 1320, 3674, 10318 






1, 0, 0, 2, 0, 0, 16, 0, 0, 192, 0, 0, 2816 (A006335) 
1, 3, 14, 67, 342, 1790, 9580, 52035, 285990 






1, 0, 3, 4, 26, 80, 387, 1596, 7518, 34656, 167310 



Table 2: The group G{S) is isomorphic to D3. These five models 
have a D-finite generating function, and the last three even have 
an algebraic generating function. 
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# 


G{S) 


S 


9(1, l;n) 
9(0, 0;n) 


References 


1 


{x,y), {yx,y), 
{yx,yx'^),{x,yx'^), 
{x,y), {xy,y), 
{xy,x'^y), {x,yx'^) 




1, 1, 3, 6, 20, 50, 175, 490, 1764, 5292 (A005558) 


Prop. El [281 [mm] 


1, 0, 1, 0, 3, 0, 14, 0, 84, 0, 594, 0, 4719 (A005700) 


2 


{x,y), {xy,y), 
{xy,x^y), {x,x^y), 
{x,y), {xy,y), 
{xy,x'^y), {x,yx^) 


/ 


1, 2, 7, 21, 78, 260, 988, 3458, 13300, 47880 (A060900) 


[MUSIBI] 


1, 0, 2, 0, 11, 0, 85, 0, 782, 0, 8004, 0, 88044 
(A135404) 



Table 3: The group G{S) is isomorphic to D^. Both models have 
a D-finite generating function, and the second one even has an 
algebraic generating function. 




Table 4. These 56 step sets are associated with an infinite group. 
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Fixed point 


Condition 




Four steps 


E 


(a, a ■2) 


a* + - 1 


X'^ + gjf'" + 31X^ + 62X^ + nx"- + 62X^ + 31X^ + 9X + 1 


H 


(«.^) 


- a'' + + 2 - 1 


X"' + 18 X'' + 125 +439X'' +897X'' + 1131 X^ +897X* + 
439 X^ + 125 X^ + 18 X + 1 





(a, a) 


a" + - 1 


X'^ - igx'^ - X" - 124X5 _^ 3jf4 _ i24X^ - X^ - 19X + 1 





(672, fe) 


6'^ + - 4 


X'^ + 9 X^ + 59 X** + 70 X^ + 59 X^ + 9 X + 1 


Five steps 


m 


(-1,1) 




2X^ + 3X + 2 


m 




+ + 2a^ - 1 


X^" + 16 X" + 106 X^ + 371 X'' + 764 X^ + 967 X^ + 764 X^* + 
371 X^ + 106 X^ + 16 X + 1 


E 


(a, a) 


a'^ — a — 1 


X*^ + 8 X^ + 28 X* + 41 X^ + 28 X^ + 8 X + 1 





(a, a^'^) 


- a - 1 


X''+2X5+6X'*+5X^+6X^+2X + 1 


m 




36® + + &2 - 1 


X" + 14X5 + g7 jf4 j^QQ ^3 + 87 + 14 X + 1 





(a, 1/a) 


+ a + 1 


X'' + 17X2+81X^ + 17X + 1 


Six steps 





(a,2 + a-3a2-3a3) 


4a^ + 1 


X^^ + 23 X^l + 283 X^" + 1861 X^ + 7461 X^ + 14225 x'' + 
18249 X'^ + 14225 X^ + 7461 X'* + 1861 X^ + 283 X^ + 23 X + 1 





/ 2+26-6^ ,\ 


fe^ + 6*' + 26^ + Sb"* + 
- 4b2 - 86 - 4 


X^'* +37X1^ + 567X1^ +4853X1^ +26197X1" +89695 X^ + 
194611 X"* + 250446 x'' + 194611 x" + 89695 X^ + 26197 x" + 
4853 X^ + 567 X^ + 37 X + 1 





{a, a) 


+ - 1 


X** + 117X'' + 3671X'' + 13396X5 + 19683X* + 13396X^ + 
3671X^ + 117X + 1 


Seven steps 





(a, 1/a) 


a'^ — a — 1 


X'' + 28 X5 + 224 X* + 345 X^ + 224 X^ + 28 X + 1 



Table 5: Proving that G{S) is infinite. 
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